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Abstract. The X-ray spectrum from a cool accretion disk heated by virialized protons is computed. The cool disk 
is either embedded in a magnetically heated accretion disk corona or partly extends into an ion supported torus 
(or ADAF). We calculate the stationary equilibrium between proton heating, electron thermal conduction and 
the radiative losses by bremsstrahlung and Compton scattering. A heated surface layer on top of the accretion 
disk is produced with temperatures between 60-90 keV above a cool layer with temperatures of 0.01 keV (AGN) 
and IkeV (galactic black hole candidates). The spectra produced by the surface layer are reminiscent of hard 
state spectra, but a bit too steep, especially for AGN's. Near the inner edge of the disk, where the optical depth 
of the disk r < 1, we find that the cool component of the disk disappears. Instead, the hot protons from the 
corona/ADAF heat the disk, on a dynamical time-scale, to temperatures of several 100 keV, limited by pair 
production. This region, here called a 'warm disk', could contribute significantly to the hard X-ray spectra and 
could be important for feeding material into an ADAF. 
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1. Introduction 

The X-ray spectra of galactic black hole candidates (BH) 
and active galactic nuclei (AGN) are remarkably similar, 
in spite of their large differences in mass and length scales. 
The similarities of the observed spectra are evident in their 
power law shape of energy index s ~ 1 in the medium 
X-ray range, and their high energy cut-off at Ec ~ 200 
keV. It is widely accepted that the power law spectra are 
produced by inverse Compton scattering of soft photons 
on hot thermal electrons (Shapiro et al. 197(:; sunyaev 



state, with the ions near their virial temperature in the 
inner region of the accretion disk. This model could in 
principle account for the observed Comptonized spectra. 
In this model the accretion flow is optically thin and ra- 
diatively inefficient. It was called an ion supported torus 



by Rees et al. (1982) 



& Titarchuk 1980; Pozdnyakov et al. 1983). Such spectra 
depend on the optical depth and the temperature of the 
Comptonizing region. In view of the similarities in the 
observed power laws these two parameters apparently vary 
only within narrow ranges for different sources. 

Another common feature evident in many of the ob- 
served spectra is the signature of a fluorescent Fe emis- 
sion line and a Compton-reflection component. These two 
components provide a strong indication that the X-ray 
production region is very close to cold matter in the cen - 
tral regions of BHs or AGNs (e.g. Ross fc Fabian 1993 ). 



The solution described by Shapiro et al. (1976 ) turned 
out to be thermally unstable. Stable accretion is obtained 
when the advection of internal energy with the flow is 
taken into account ( Narayan fc Yi 1994 , 1995a| ,[b|); these 
models are called advection dominated accretion flows 
(ADAF). Most of the viscously dissipated energy is ad- 
vected radially with the flow. The protons are near their 
virial temperature, whereas the electrons are much cooler 
due to their strong interaction with the radiation field and 
the low rate at which they can exchange energy with the 
protons via Coulomb interactions. The flow may in prin- 
ciple coexist with an optically thick, cold accretion disk 
such that the cold disk partly extends into the hot torus 



(Esin et al. 1997) 



The classical accretion disk model ( Shakura & Sunyaev 



and 



Another possibility for the simultaneous presence of 
cold and hot matter are accretion disk corona (ADC) 



1973) can not explain the simultaneous presence of cold 

h( tt matter near a compact object. Another mo de 1991 



Nakamura & Osaki 1993; 


Haardt & Maraschi 


Svensson & Zdziarski 199^; 
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Dove et al. 1997). 



of acc retion was subsequently proposed by Shapiro et al 
1976| ). They showed that there is an optically thin hot 



disk solution where the plasma is in a two-temperature 
Send offprint requests to: bed@mpa-garching.mpg.de 



In these models the cold accretion disk is embedded in a 
hot corona in a plane-parallel slab conflguration. A large 
fraction of the gravitational energy is assumed to be re- 
leased in the corona by magnetic fields, although the de- 
tails of the coronal heating process are still unclear. The 
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protons in the corona are hot, whereas the disk beneath 
is cool {kT ~ 1 keV), optically thick, geometrically thin 
and gas pressure supported. 

In both the ADAF and the ADC models there is an 
energy coupling between the cold disk material and the 
hot tenuous plasma above. Traditionally this interaction 
was seen in terms of an exchange of radiation ( Haardt & 
Marasdhi 19911 



1993 



Haardt et al. 1994, 1997). Here we 



are interested in the energetic coupling via the exchange 
of matter. We investigate the penetration of hot protons 
from the ADC/ ADAF into the cool disk, i.e. 'ion illumina- 
tion'. The hot protons are stopped in the cold disk due to 
Coulomb interactions with the ambient electrons, rather 
than by binary interactions with target ions. 

Proton illumination is not restricted to accretion disks. 
It was applied early in X-ray astronomy to explain spectra 
of galactic X-ray sources as neutron star surfaces heated 



by radially infalling ions ( 


Zel'dovich & Shakura 1969; 


Alme & Wilson 1973 


; TuroUa et al. 1994 


). Recently 


Deufe^ 



et al. (;!001) revisited the proton illumination of a neutron 
star in the context of accretion from an ADAF. 

The importance of proton illumination for accretion 
disks was su ggested by ^pruit (1997 ) and Spruit fc 



Haardt (2000). Detailed Monte Carlo calculations of the 
Comptonization were presented by Deufel & Spruit (2000) 
[henceforth paper I]. For the present work we have im- 
proved the treatment of the radiative processes by solving 
the radiative transfer equation including thermal emission 
due to bremsstrahlung. Pair production is also taken into 
account. We further allow for energy redistribution due to 
thermal electron conduction within the disk. The inter- 
action between the impinging protons and the accretion 
disk is computed time dependently in a one-dimensional, 
plane-parallel approximation. The density distribution is 
found from hydrostatic equilibrium including the force ex- 
erted by the decelerating protons. 

The cool disk acts as an effective thermalizer for ener- 
getic (Comptonized) photons as long as its optical depth 
is large enough. At lower optical depth, such as occurs 
near the inner edge of the disk, photon production is less 
efficient. Therefore we expect different solutions for disks 
of low optical depth. This is borne out by the results re- 
ported in Section |5[ 

2. A cool disk extending inside an ion supported 
flow 

2.1. Optically thick case 

Most of the physics of the interaction between the hot sup- 
ported flow and the underlying cool disk is rather indepen- 
dent of the structure of the cool disk, since it acts mainly 
as a thermalizing surface. The height above the midplane 
where the hot ion-supported flow and the cold disk inter- 
act is important, however, since this determines the grav- 
itational acceleration, and hence the density within the 
interaction layer. Processes such as photon production by 
bremsstrahlung depend on the electron density explicitly. 



To obtain an estimate for the thickness of the cool disk, the 
amount of energy released in the cool disk must be speci- 
fied. For this purpose we adopt an a-disk prescription as 
described in Svensson and Zdziarski (1994). In this model, 
a fixed fraction / of the gravitational energy release takes 
place above the disk, and a fraction 1 — / in the cool disk. 
For details we refer also to paper I. 

In order to simplify our treatment we assume that the 
vertical disk structure can be approximated at a given 
radius by a plane-parallel, one-dimensional geometry, so 
that the vertical z-coordinate is the only dimension rele- 
vant to our problem. We find the vertical density distri- 
bution from the equation of hydrostatic equilibrium. 



dP 



K OT 



(1) 



1/2 

where — [GAI/R^^ is the local Kepler angular ve- 
locity, K = Kos = 0.40 cm^g~^ is the electron scattering 
opacity and r is the electron scattering optical depth. Note 
that we use r here as a coordinate, replacing the verti- 
cal height z. The first term on the right represents the 
usual weight of the disk; the second term takes into ac- 
count the momentum exerted by the deceleration of the 
incident protons as a function of optical depth. The force 
dPp/dr is evaluated by recording the change of velocity 
Avp of the incoming protons as a function of optical depth 
(see Sect. 3T). The contribution of the radiative pressure 
is neglected as we study cases well below the Eddington 
luminosity. 

Together with the equation for the vertical coordinate 



dz 
rf7 



~k T{t) 
firup K P{t) 



(2) 



we integrate Eq. (|i|) via a fourth order Runge-Kutta 
method. We set /Lt = ^ for the idealized case of an ion- 
ized hydrogen atmosphere. 

The initial value for the pressure at the top of the 
layer is set to a small fraction of the coronal pressure, 
-ftop(''' = 0) = lO^^ripfcTvir. The geometrical height ztop 
of this upper boundary layer is not known in advance. We 
solve the pressure profile by starting with an initial guess 
for the height of the disk above the mid-plane, Ztop- We 
integrate Eq. (|l|) to the maximum optical depth of our 
model atmosphere Tbot at the vertical height Zbot- Usually 
we set Tbot = 3. At this Thomson optical depth the pro- 
tons have already lost more than 0.99 of their energy. At 
''bot pressure balance between the simulated layer and the 
underlying cool disk requires -P(Tbot) = P{zhot)- We iter- 
ate until we find the right value for which this pressure 
condition at Zhot is fulfilled. Thus we have matched the 
simulated layer to the underlying cool disk. The temper- 
ature of the cool disk and the pressure at the mid-plane 
are determined as in paper I. 
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2.2. Moderate optical depth 

Near the inner edge of a cool Shakura-Sunyaev disk at the 
radial distance Ri from the compact object, the surface 
density drops to low values. In the limit H/R ^ 1 (where 
H is the disk pressure scale height), the surface density 
varies as [1 — Close to the inner edge of an 

accretion disk, the assumption of a large optical depth is 
therefore not valid any more. A certain minimum optical 
depth is necessary, however, for efficient cooling of the disk 
by soft photons to be possible. In section || we explore 
cases where the total optical depth of the layer is of order 
unity or less. The procedure to calculate the hydrostatic 
balance in such cases is almost the same as in the cases of a 
surface layer on top of an optically thick disk. But instead 
of matching the computed layer to an underlying cool disk, 
we put the bottom of the simulated layer is at 2;(rbot) — 0, 
the mid-plane of the disk. We are interested in simulating 
the complete disk, i.e. from the upper surface to the lower 
surface. We do this by introducing an artificial 'mirror' 
placed at the mid-plane. All relevant physical quantities 
of the model are reflected at this mirror. Thus we have 
to compute one half space only instead of simulating the 
whole vertical disk structure. 



3. Heating and cooling processes of the model 

3.1. Proton illumination of the accretion disk 

3.1.1. Coulomb interactions in an ionized plasma 

At the typical energies of the protons incident on the cool 
disk, the energy loss is mostly by long-range Coulomb in- 
teractions with the electrons in the disk (small-angle scat- 
tering on the large number of electrons in a Debye sphere). 
This is opposite to the case of protons with a tempera- 
ture near that of the plasma in which they move. In the 
latter case, the equilibration among the protons is faster 
than between electrons and protons, by a factor of order 
(mp/mc)^^^. To see how this apparent contradiction is re- 
solved, consider the basic result for the energy loss of a 
charged particle moving in a fully ionized, charge-neutral 
plasma. This was derived by Spitzer (1962) (making use of 
Chandrasekhar's (1942) earlier results on dynamical fric- 
tion). Introduce as a measure of distance in the plasma 
the Thomson optical depth r, i.e. dr = (TttkU, where 
(7t = 87re*/ (3mgC^) is the Thomson cross section, n the 
electron density and I the distance. The rate of change 
of energy E — TOpU^/2 of a proton moving with velocity 
V in a field of particles with charge e and mass nif (the 
'field particles' in Spitzer's nomenclature) is then given by 
Spitzer's Eq. 5-15. In our notation, this can be written in 
terms of the energy loss length Tf for interaction with the 
field particles f , 



^f"' = = 31n A(^)2(^)4(l + ^)F(xf), (3) 

h dr nip v m{ 



where In A is the Coulomb logarithm (which is determined 
by the size of the Debye sphere). Here 



F{x) — ip{x) — xip'{x), 



(4) 



where ip is the error function and ip' its derivative, and 
Xf = v[mf/{2kT)]^/^ is (up to a numerical factor) the ratio 
of the incident proton's velocity to the thermal velocity of 
the field particles. The limiting forms of F are 



F{x)^x^/3 (x^O), 



1 (a; — > oo). 



(5) 



We can evaluate (|^) under the assumption that the field 
particles that are most relevant for the energy loss are the 
protons or the electrons, respectively, and compare the loss 
lengths. If the incident proton has velocity comparable 
with the thermal velocity of the field protons, we have 
Xp ^ 1, and Xc — {mc/mpY^'^Xp <^ 1. F{xp) is then of 
order unity and F{xc) ~ x^/'^- Setting f = e respectively 
f — p in (||) and taking the ratio, we have 



— ~ / "^P \l/2 
Td Too 



(6) 



The loss length for interaction with the electrons is thus 
much longer than for interaction with the protons, and 
the interaction with electrons can be neglected. This well 
known result is the relevant limit for the relaxation of a 
proton distribution in a plasma that is not too far from 
its thermal equilibrium. 

For incoming protons of high energy, however, the re- 
sult is different because Xc is not sufficiently small any 
more. In the high-t; limit, F{xc) ~ F{xp) — 1, and one 
has 



2— = 10" 

TOn 



(7) 



In this limit, the energy loss is thus predominantly to the 
electrons. A related case is that of the ionization losses of 
fast particles in neutral matter (for references see Ryter et 
al. 1970). The case of an ionized plasma is simpler, since 
the electrons are not bound in atoms. The change from 
proton-dominated loss to electron-dominated loss takes 
place at an intermediate velocity Vc, at which F(xp) « 1 



but Xc still small, so that F(xc) 
Tp then yields 



^ Too 
TOn 



x'^/3. Equating To and 



(8) 



Ec 
kT 



12. 



(9) 



For the electron temperatures we encounter in our models, 
T ~ 100 keV, energy loss to the field protons can thus be 
neglected for incoming protons with energy E ^ 1 MeV. 
This is the case in all calculations presented here. 
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3.1.2. Corrections at high and low energies 

Spitzer's treatment is non-relativistic, while virialized 
ADAF protons near the hole can reach sub-relativistic 
temperatures. A fully relativistic treatment of the 
Coulomb interactions in a plasma has been given by 
Stepney and Guilbert (1983). We have compared the clas- 
sical treatment according to Spitzer's theory with this rel- 
ativistic result in Deufel et al. (2001), and found it to 
be accurate to better than 5% for proton temperatures 
< 100 MeV. The classical approximation in Spitzer's anal- 
ysis therefore does not introduce a significant error for the 
problem considered here. 

For high energies, the Coulomb energy loss becomes so 
small that loss by direct nuclear collisions becomes com- 
petitive. This happens (cf. Stepney and Guilbert 1983) at 
E ^ 300 McV, an energy that is not reached by virialized 
protons except in the tail of their distribution. We ignore 
these direct nuclear collisions. Note, however, that a grad- 
ual nuclear processing by such collisions can be important 
(Aharonian and Sunyaev, 1984), in particular for the pro- 
duction of the Lithium. The Lithium overabundances seen 
in the companions of LMXB (Martin at al., 1994a), may 
in fact be a characteristic signature of the interaction of 
an ADAF and a disk described here (Martin et al. 1994b, 
Spruit 1997). 

As the protons slow down, they eventually equilibrate 
with the field protons. This last part of the process is not 
accurately described by the energy loss formula (^). In 
addition to the simple energy loss of a particle moving on 
a straight path through the plasma, one has to take into 
account the random drift in direction and energy resulting 
from the interaction with the fluctuating electric field in 
the plasma. This drift can be ignored to first order (end 
of section 5.2 in Spitzer 1962), but takes over in the final 
process of equilibration with the plasma. This last phase 
involves negligible energy transfer compared with the ini- 
tial energy of the protons in the present calculations, and 
can be ignored here. 

3.1.3. Charge balance 

The protons penetrating into the disk imply a current that 
has to be balanced by a 'return current'. As in all such sit- 
uations, this return current results from the electric field 
that builds up due to the proton current. This field drives 
a flow of electrons from the ADAF to the disk which main- 
tains the charge balance. Since the electron density in the 
disk is high, the return current does not involve a high 
field strength. 

3.1.4. Calculations 

The numerical method used to calculate the electron heat- 
ing rates of the accretion disk due to proton illumina- 
tion is described in detail in paper I; here we give a short 
overview. 



We follow the evolution of an initially Maxwellian dis- 
tribution of protons, which we place above the cool ac- 
cretion disk, as the protons penetrate through the disk 
atmosphere. For the temperature of the incident protons 
in our model we take the local virial temperature 



n — J- v 



GMnip 
3kBR 



156 



MeV, 



(10) 



where G is the gravitational constant, M the mass of the 
black hole, mp the proton mass, R the radius from the 
central object and is the Boltzmann constant, r is the 
dimensionless radius, r = R/Rs, where Rg — 2GM/c^ is 
the Schwarzschild radius. 

The energy loss of the incoming protons is computed 
from (^), with f = e. Denote by r the Thomson optical 
depth measured vertically into the layer, from its top. The 
rate of change of energy with depth for a proton at depth 
T moving at an angle 9 with respect to the vertical is then 



1 AE 31nATOe c 4 
E dr cost' m„ v 



with 



X = Xc = v(——) 



1/2 



(11) 



(12) 



while T(t) is the depth-dependent temperature of the 
layer, and v the depth and angle-dependent velocity 
of the incident proton. The angular distribution of 
the protons is discretized in 50 equidistant points. 

The Debye-length in the plasma, measured in Thomson 
optical depths, is quite small so that temperature varia- 
tions over a Debye-length can be ignored and the local 
energy loss rate at any depth is adequately described by 

The energy fiux of the incident fiux qp is proportional 
to the density in the ADAF, which depends on its accre- 
tion rate and viscosity. Instead of using detailed models 
we parameterize these dependences by scaling the energy 
flux of the incident protons with the local energy dissipa- 
tion rate in the ADAF. Using the thin-disk expression for 
this dissipation rate, we set 



,„^/<3(«) = /-5||M.„«,, 



(13) 



where / is an adjustable numerical factor of order unity 
and J{R) = 1 - (Ri/R)^/'^. This scaling of the proton flux 
is also useful if one has a coronal model in mind instead 
of an ADAF. The factor / then is interpreted as the frac- 
tion of the viscous energy release that is dissipated in the 
corona. 

We follow the protons from the corona numerically 
through the disk atmosphere and record the energy loss as 
a function of optical depth r according to Eq. (|ll|). This 
yields the local time dependent heating rate Ap(t). 
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3.2. Radiative transfer 

We treat the radiative transfer by solving the radiative 
transfer equation 



dl 



xiX2d\og{y'WiX2)d\og{xi / X2) , since the function G{xiX2) 
is independent of xxjx^- 

The pair production rate due to electron-electron col- 
lisions is e.g. given in ^vensson "(198^ as 



dz 



s 



(14) 



(n++n_)2cr2Pee(e) 



where 1^ ^ is the intensity as a function of frequency v where we have introduced Q = kTc/rricC^ and 
and photon angle /i — cos(0), is the bremsstrahlung 
absorption coefhcient, the Compton extinction coef- 
ficient, the bremsstrahlung emissivity and jj^ ^ the 



P., 



^(2 In 20)3 
2 ■'^ 277r^ ' 



(19) 



(20) 



Compton emissivity. The method fo r solving the radia - 
tive transfer equation is described in Deufel et al. (2001 ) . 
From the solution of the radiative transfer equation we 
obtain Ittie radiative cooling rates Ai.ad('''j- 



Here a/ is the fine structure constant. This production 
rate is only included for > 0.6. 



(1982) as 



3.3. Pair processes 

As discussed below, the cases of moderate optical depth 
can become rather hot so that pair processes have to 
be included. We include photon-photon pair production 
(77 ~* e~^e~) as well as pair production due to ee collisions 
(ee — > eee+e") in a steady state, i.e. the pair production 
rate equals the pair annihilation rate. 

The photon-photon pair production rate can be ex- 
pressed as a quadruple integral over the dimensionless 
photon energies xi = hvi/rrieC^ , X2 — hv2/meC^ and pho- 
ton angles /xi, /i2, respectively (cf. Zane et al. 1995): 



The pair annihilation rate is also given by Svensson 

(21) 



where 

TT 

m = 



1 + 2612/ 111(2776! + 1.3) 



(22) 



7] ~ 0.56146. Now we can set up the pair balance equation 
in the stationary case. 



R 



.^^ + c ri [{2z + 1)^^(61) - z(l + z)A{e)] = 0. (23) 



^-77 = 7^ / dxidx2duidu,2 ^^''^l, tj^''^^ F(x+, X-), (15) 



We have used the charge neutrality condition ri_ = n++ n 
and z = n+ /n is the ratio of the pair number density to the 
proton number density n. Eq. ( p^ ) is a quadratic equation 
for the pair fraction z. The positive root of this equation 



where F{x-i-,x^) is defined as 

a{x)x^dx 



F{x+,x.)= I 



z{r) 



^ {x\ — a;2)(x2 — x'L) 



(16) 



(24) 



with x\. = a;ia;2[l — cos(0i ± 6l2)]/2 ( IStcpney fc Guilberl| 
1983). The pair cross-section cr{x) is (Jauch & Rohrlich 



197e) 



The additional electrons from pair production serve as 
further scattering partners for the Coulomb collisions and 
the Compton scattering. Effects due to electron-positron 
bremsstrahlung are not included nor is a treatment of the 
annihilation line. 



<j{x) - -{1-13^) 



(3-/3') log 



1 + /3 
1-/3 



2/?(2-/32) 



,(17) 3.4. Energy balance from heating and cooling 



with /? = y/x^ — 1/x. 

In our models we have simplified the above integral 
by replacing the intensities by their corresponding mean- 
intensities: J VI and /^2,i'2 ~^ f'o^ ^iiy given 
value of y = X1X2 we can now evaluate the function 
G{y) = J F{x+,X-)d^iid^2, which is independent of the 
radiation field, and which can be tabulated and stored 
beforehand. The pair production rate now reduces to 



« J G{xiX2)JuiJv2dxidx2 



(18) 



Thus the quadruple integral is reduced to a double inte- 
gral, which is much faster to evaluate. An accurate eval- 
uation of this integral can be made by writing dxidx2 = 



We begin our computation with an isothermal tempera- 
ture profile in hydrostatic equilibrium according to Sect. 
||. The initial temperature throughout the layer is Td (see 
paper I) for the thick disks and we set Tc = 1 keV for 
disks with moderate optical depths. The solutions do not 
depend on the initial temperature profile. First we calcu- 
late the heating rates A+(t) from the Coulomb interac- 
tions and the cooling rates A~^(r) due to the radiative 
processes bremsstrahlung and multiple Compton scatter- 
ing. Additionally the energy redistribution due to electron 
thermal conductiv ity is included using S pitzer's classi- 
cal value, as in [Deufel et al. (200lj ). This process 



adds the contribution Acond(T) to the energy bal- 
ance. The validity of classical electron conduction 
was checked by evaluating the electron mean free 
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path A, which turns out to be of the order 10^ cm or 
less in the black hole candidate cases shown in Fig. 
|l|. This is much smaller than the temperature gra- 
dient length (more than 10^ cm), so the condition 
for validity of the electron formula used is amply 
satisfied. The same holds for the AGN cases. 

The equilibrium is computed following the time evo- 
lution of the layer until a balance between the heating 
and cooling processes is obtained. In the optically thick 
models, the time scales for approaching equilibrium turn 
out to be a sharp function of depth in the model. To deal 
with this, an adaptive time stepping process is used in 
which the time step depends on both time and depth in 
the model. Stability of this process was obtained by scal- 
ing this step with the shortest of the energy exchange time 
scales associated with the contributing heating and cool- 
ing processes. With this procedure the time evolution of 
the model is not realistic, but the final equilibrium ob- 
tained is. For the cases with moderate optical depth we 
use a depth-independent time step. In these cases, the time 
evolution of the model is also realistic. 

The total change of enthalpy per time step as function 
of optical depth can then be expressed by 



Aw{t) 
At 



pcp- 



AT{t) 



At 



AjM+A-^(T)+Aeond(T),(25) 



where Cp is the specific heat at constant pressure. 

Now the change of temperature as a function of optical 

depth CPiri bp palrnlatpH AAA'tVi tViP np-^sr tpmppratnrp prn- 



The upper right panel of Fig. ^ shows the heating and 
cooling rates. In the narrow transition zone (Ar « 0.1 ) 
between the hot top and the cool bottom electron con- 
ductivity is not negligible. Elsewhere conductivity plays 
no role for the energy balance. 

The geometrical depth of the heated layer is ztot ~ 
6x 10^ cm. This is small compared to the distance from the 
compact object {R = 2.3 x 10^ cm). Thus a cool disk with 
a proton heated skin on top of it can still be considered a 
"thin' disk. 

Fig. ^ shows the dependence of the emergent spectrum 
and the temperature profile on the distance from the com- 
pact object in the galactic BH case. We obtain a result 
comparable to paper I, i.e. with increasing distance the 
temperature of the hot part increases somewhat, whereas 
the optical depth of the heated layer decreases. A discus- 
sion of these effects can be found in paper I. 

For distances R > 50i?s from the galactic BH we could 
not find heated solutions. With increasing distance from 
the central object the proton temperature (cx R~^) as well 
as the proton energy flux (cx R^^) decrease. In these con- 
ditions we find solutions in which the energy supplied by 
the incident protons is thcrmalizcd directly in the cool disk 
by bremsstrahlung, without forming a hot surface layer. It 
is likely that a very thin (r ^ 0.01) hot layer still forms in 
these conditions, but we are unable to resolve it with the 
present method. This is indicated by the existence of thin 
hot atmosphere solutions at low incident energy fiux in the 
case of protons heating a neutron star surface (Zampier: 



file tViplViyrlrngtatip titriirtnrp ig iipHatpH arrnrHing tn Sppt 

^. We follow the simulation until the Coulomb heating is 
balanced by the radiative cooling and energy redistribu- 
tion due to electron conductivity. The energy balance of 
all our computations is better than 10%. 

4. Results of the model computations 

4.1. Proton illumination of cool disks in BH and AGN 

For the thick disks we use the same model parameters as 
in paper I, i.e. the accretion luminosity L — O.lLEdd, the 
viscosity parameter a = 0.1 and the fraction of the energy 
released in the virialized atmosphere / = 0.95. The mass 
of the galactic BH is A/bh — 8Mq and the mass of the 
AGN is A/agn = 8 X IO'^Mq. 

Fig. |l| shows the solution at r = 10 for a proton energy 



et al. 199£; Deufel et al. 2001). It probably also occurs in 



flux of Fp = 8.1 X 10 erg cm 



T„ = 15.6 MeV. The 



cool disk is separated by a sharp temperature front from 
a hot part at an optical depth of order unity. This step 
like temperature profile is much more pronounced com- 
pared to the results of paper I. The temperature of the 
hot part is Tc ~ 60 keV. The hot layer acts as an effective 
Comptonizing region and a hard spectrum is emitted. The 
power law index of the BH spectra in EF{E) is s 0.2. At 
the soft end of the spectrum there is an excess of soft pho- 
tons with respect to a blackbody with equivalent energy 
flux. This exces s is due to a 'rever se photosphere effect', 
as explained in Deufel et al. (2001). 



the present case, but since these thin layers at large dis- 
tances from the hole do not contribute much to the overall 
spectrum, we have not pursued this further. 

Fig. H shows the dependence of the emergent spectrum 
and the temperature profile on the distance from the com- 
pact object in the AGN case. We obtain slightly different 
solutions compared to our results in paper I. The optical 
depth of the heated layer is smaller (r < 0.5), the tem- 
perature of the hot part is Tq k, 70-90 keV. Again a steep 
temperature jump separates it from the cool disk under- 
neath. Due to the smaller optical depth of the heated layer 
Comptonization is not as efficient. A significant fraction of 
the incident ions passes through the hot layer and thermal- 
izes in the cool disk underneath. This increases the soft 
photon flux relative to the Comptonized flux, and pro- 
duces a steeper spectrum. The index of the AGN spectra 
in EF{E) is s « 0.6. Our model does not produce heated 
solutions at distances R > 20i?s from the compact ob- 
ject, though heated atmospheres might exist at very small 
optical depths (see note above). 

5. Proton illumination of disks with moderate 
optical depth 

The temporal evolution of an initially cool (To — 1 keV) , 
thin disk with an optical depth — 1-5 (measured from 
the surface to the mid-plane) at r = 10 for the galactic 
black hole candidate (BHC) case is shown in Fig. ||. The 
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Fig. 1. X-ray spectra of proton- illumi- 
nated optically thick cool disks. Upper 
left panel: emergent model spectrum 
(solid line) and blackbody with tem- 
perature according to proton energy 
flux (dotted-dashed line); upper right 
panel: combined heating rates from pro- 
ton heating and electron thermal con- 
ductivity (solid line) and radiative cool- 
ing rates due to Comptonization and 
bremsstrahlung (squares) - the dotted 
line shows the rates from electron con- 
ductivity alone; lower panels show from 
left to right and from top to bottom 
the electron temperature To, pressure 
P, electron density Wc and the geomet- 
rical depth of the layer z/ R, for a solu- 



tion with F„ 
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Fig. 2. Dependence of the emergent model spectra and the temperature profiles on the distance from a galactic BH 
with M = SMq. The surface temperature of the of the disk increases slightly with distance, whereas the optical depth 
of the hot surface layer decreases. 



energy flux from the virializcd protons is Fp — 8.1 x 10^^ 
erg cm~^ s~^. In the beginning the temperature of the top 
layers increases due to the impinging hot protons, whereas 
the mid-plane region cools due to bremsstrahlung [stage 0- 
2, Fig.|l. As the top layers are heated, the stopping power 



of the plasma decreases and the protons penetrate deeper 
into the disk. Eventually hot protons reach the mid-plane 
region and proton heating overcomes bremsstrahlung cool- 
ing even there [3] . At this stage no effective cooling mech- 
anism is present and the temperature continues to rise ev- 
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Fig. 3. Dependence of the emergent model spectra and the temperature profiles on the distance from the hole in the 
AGN case with Af = 8 x IO^Mq. The optical depth of the hot part is much smaller than in the galactic BH cases. 
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Fig. 4. A 'warm disk', of optical depth 
1.5, externally heated by ion illumina- 
tion. Upper left panel: emergent model 
spectrum; upper right panel: heating 
rates from proton heating (solid line) 
and combined radiative cooling rates 
from Comptonization, bremsstrahlung 
and pair processes (squares); lower left 
panel: electron temperature T^; lower 
right panel: proton number density 
(solid line) of the layer and pair number 
density (dashed line). 



erywhere in the disk [4-5]. At temperatures kT > 200 keV 
pair production becomes more and more important. The 
extra electrons serve as additional scattering partners for 
the Coulomb collisions with the protons and the pho- 
tons of the radiation field. Thus pair production limits 
the maximum attainable temperatures, and the disk ad- 
justs to a new equilibrium state [6] at a temperatures 
Te ~ 400 keV. The dashed line in Fig. || shows the temper- 
ature profile after one dynamical time-scale at that radius, 
td = l/^K ^ 3.5 X lO-^s. 

Fig. H shows an overview of that solution. The pair 
number density z at the top of slab reaches z « 50% of the 
proton number density and drops off at the mid-plane to 
z « 1%. The spectrum of such a thin proton heated disk 
peaks at w 1000 keV. We refer to those disks as warm 



disks . They are still considerably cooler than the local 
virial temperature. 

Fig. ^ shows the dependence of the equilibrium tem- 
perature profile on the initial optical depth of the layer 
and the distance from the central object. With increasing 
distance the proton energy flux as well as the proton pene- 
tration depth decreases. A cool disk can therefore be trans- 
ferred into the warm state only within a certain distance 
from the BH. If the hot protons do not reach the mid-plane 
anymore, a cool interior can be maintained which looses 
its energy very efficiently via bremsstrahlung, as the hot 
layers above are optically thin. Our model computations 
show that at r = 15 a thin disk can maintain a cool inte- 
rior for Ti/2 = 1-5 whereas at r = 10 a disk with Ti^2 ~ 1-5 
switches into the warm state. 
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Fig. 5. Evolution from a cool disk to a warm disk by ex- 
posure to virialized protons at r = 10. Initial temperature 
Tc = IkeV The numbers at the lines denote the sequential 
stages, referred to in the text. The final stage [6] shows an 
equilibrium at Tc w 400 keV. The dashed line shows the 
temperature profile after one dynamical time-scale. 




0.01 



Fig. 6. Equilibrium temperature profiles of warm disks 
exposed to virialized protons at r = 10 (solid line) and 
r = 15 (dashed line) with different initial optical depths 
Ti/2 = 2, 1.5, 1. for r = 10 and = 1-5, 1 for r = 15. 



The temperature of the warm state also depends on 
both the distance and the overall optical depth of the 
layer. For T1/2 = 1 and r = 10 our model predicts a tem- 
perature of w 1 MeV. At such temperatures our classical 
proton heating formalism starts deviating from the correct 
relativistic expression. The classical treatment underesti- 
mates the proton-electro n heating rates at high tempera- 

But further pair processes and 



tures (Deufel et al. 2001 



radiative cooling terms should also be included (see be- 
low), which again limits the maximum temperatures. 

The transition from a cool disk to the warm state also 
takes place in AGNs. For the above AGN parameters (see 
Sect. O) we find the transition to occur for ^ 0.2 
at r < 10. The temperatures of the warm state are in the 



MeV range. At such temperatures {6 > 1) our treatment 
of the pair processes and the radiation field needs more 
scrutiny. Further pair production processes should be in- 
cluded (76 ee^e~ , ep — > epe'^e") as well as additional 
radiative cooling terms (bremsstrahlung from e~^e~' and 
e^e^ collisions). These will further limit the maximum 
attainable temperatures. Therefore we do not think that 
temperatures of several MeV found in our simulation of 
the warm state in AGNs are realistic. The important re- 
sult of our investigation is that the transition is not only 
restricted to the BHCs but also takes place in AGNs. But 
compared to the BHCs the transition in the AGNs oc- 
curs in a more narrow zone around the BH (in terms of 
Schwarzschild radii) and the vertical extend (in terms of 
Thomson optical depths) is smaller. 

6. Discussion and conclusions 

We have investigated a model for the production of a hard 
spectral component {hv ^ IkeV) from a cool disk around 
a compact object. The model assumes that a standard 
cool accretion disk is surrounded by a hot atmosphere. 
This atmosphere can either be due to an ion supported 
accretion flow (ADAF) which partially overlaps in radial 
extent with the cool disk or to a hot corona which is e.g. 
magnetically heated (e.g. |di Matteo et al. 1999 ). For the 
actual calculations presented here, we have used parame- 
ters representative for ADAFs. Fast protons from this at- 
mosphere interact with the cool disk material and produce 
a warm (^ 80 keV) surface layer, where soft photons from 
the cool disk are Comptonized. We have considered two 
cases in detail: a disk around an SAf© black hole (galactic 
BHC case), and an AGN case with an 8 • 10^ Mq black 
hole. 

We flnd that the response of the cool disk to proton 
illumination depends on its optical thickness. At sufficient 
thickness (r > 1.5) the cool disk is able to effectively re- 
process incoming energy fiux into soft photons. Here, a 
two-layer structure develops, consisting of a warm surface 
layer of depth r ~ 1 on top of a cool disk. At lower total 
optical depth, the soft photon emissivity is insufficient to 
cool the disk and the disk heats up, by the proton energy 
fiux, to an approximately isothermal layer of intermediate 
(a few 100 keV) temperatures. We call these cases warm 
disks: their temperature is higher than an accretion disk 
radiating as a black-body, but still much lower than the 
virial temperature. 

In the optically thick case, the energy balance between 
the warm surface layer and the cool disk is like that in 
the Haardt-Maraschi (1991, 1993) model of a corona over 
a cool disk. For geometrical reasons about half of the 
Comptonized photons is reprocessed by the cool disk, the 
other escapes. The Compton y— parameter adjusts such 
that the amplification factor matches this requirement. 
This balance regulates the temperature of the layer, while 
its optical depth is determined by the penetration depth 
of the protons. The balance is strongly self-regulating: the 
temperature and optical depth of the warm layer depend 
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little (within limits) on temperature and energy flux of the 
incident protons, so that also the shape of the emergent X- 
ray spectrum varies little (see Haardt and Maraschi 1993, 
and paper I for a more detailed discussion). 

We find that the existence of these proton-heated sur- 
face layers is limited to the innermost regions r < 50 of the 
disk. At larger distances from the hole, the optical depth 
of the warm layer drops rapidly, and its contribution to 
the spectrum becomes small. 

For our BHC case the temperature and optical depth 
of the warm layer are in the right range to produce spectra 
with the main features of black hole binaries in their hard 
states. In the AGN case, the Comptonization is weaker, 
and the hard spectral component somewhat steeper (pho- 
ton index around 2.6) than in the BHC case (index around 
2.2). In both the BHC and AGN cases the spectra cut oS 
around 200 keV. The maximum distance from the hole 
where a significant warm layer is formed is somewhat 
smaller in the AGN case (i?max/-Rs ~ 20) than in the 
BHC case {Rma^/Rs ~ 50). This is due to the lower en- 
ergy flux per unit of surface area, and the reduced pene- 
tration depths in the cooler AGN disks (see also Deufel et 
al. 2001). 

Near the inner edge of a cool accretion disk (at the dis- 
tance i?i from the hole), its optical depth becomes small 
(in the thin disk limit, it varies by the well known fac- 
tor / = [1 - (Ri/R)^/^] ). Though the optical depth of a 
cool accretion disk is generally substantial at the accretion 
rates inferred for BHC and AGN, there is always a region 
of low optical depth near the inner edge. This is also where 
the disk would be exposed to the largest proton energy fiux 
from an ADAF or corona. At optical depths ti/2 < 1.5 we 
find that an initially cool disk heats up to a new equilib- 
rium with nearly uniform temperatures of a few 100 keV 
(Fig. 5). The transition takes place in a time of the order 
of the dynamical time scale l/fl. Initially, the protons are 
stopped in the cool layers before reaching the midplane. 
Due to the thermal expansion of the plasma the electron 
densities drop, and as a consequence the proton heating 
finally overcomes the cooling by bremsstrahlung (~ n^). 
A new equilibrium is found when pair production sets in. 
The pairs add to the optical depth of the layer and increase 
the radiative loss of the plasma (by Comptonization) until 
it matches the input by proton heating. 

Our solutions are local equilibria for conditions at a 
given distance from the hole. In disks in which the optical 
depth is low (< 1) near the inner edge but at least a few at 
larger distance, proton heating would produce our warm 
disk structure of approximately uniform temperature in a 
(possibly narrow) zone near the inner edge. Further out 
it would instead produce a warm surface layer overlying 
the cool disk. Both regions would contribute to the overall 
X-ray spectrum. 

Since the energy flux by proton heating is highest near 
the inner edge, the warm disk component at that location 
could contribute significantly to the overall spectrum even 
if its radial extent is limited. In our solutions, it would add 
a very flat component (photon index around 1, up to about 



IMeV, Fig. 4), so the overall spectral index could be sig- 
nificantly less than 2, as observed in some of the hardest 
spectra. The warm surface layers over the cool disk do not 
have such flat spectra (Fig. 2), and alone can not explain 
the hardest observed spectra (with indices around 1.5). 
The combination has the potential to explain the entire 
range of observed hard spectra, but for a quantitative re- 
sult the radial structure of a proton-illuminated disk needs 
to be treated in more detail. 

In this context we note that our warm disk solutions 
ignore a potentially important cooling mechanism, namely 
Compton cooling by soft photons from a cool disk. Though 
there are no cool regions internally in the warm disk, the 
nearby cool disk at larger radii produces a soft photon fiux. 
A fraction of this flux could travel radially, illuminating 
the warm disk region. The fiux of such soft photons is 
probably small, since the disk is thin at all radii and the 
radial optical depths therefore large. For quantitative as- 
sessment a two-dimensional radiative transfer model has 
to be developed, which is outside the scope of the present 
treatment. 

The warm disk solution has a combination of temper- 
ature and optical depth just in the region where a nor- 
mal, internally heated disk cannot exist. In such disks, the 
(T, r) combination characterizing our warm disks would 
lie on the SLE (Shapiro-Lightman-Eardley) branch of ac- 
cretion fiows, which is thermal/viscously unstable. The 
combination is quite stable in our warm disk solution pre- 
cisely because it is not heated internally but externally 
by the incident protons. In this context it is interesting 
that model fits of low/hard states in terms of stationary 
ADAF models sometimes yield parameter combinations 
on the unstable SLE branch (for example in XTE 1118 
-1-480 as analyzed by Frontera et al. 2001). We interpret 
this as a strong indication for the existence of our warm 
disk solutions. 

The physics of the warm disk solutions presented here 
has further interesting consequences. In anticipation of re- 
sults to be given in more detail in a future paper, we note 
here that at the temperatures and densities in our warm 
disks, the time scale for establishing thermal equilibrium 
between the electrons and the protons in the disk (not the 
illuminating protons) is not short compared to the dy- 
namical time scale any more. Any mechanism which now 
provides energy to the disk protons faster than they can 
exchange energy with the disk electrons would give rise 
to a heating of the disk protons and, depending on the 
temperature and density dependence of such a process, 
might lead to runaway heating of the protons. Two energy 
channels immediately come to mind: proton-proton inter- 
actions, i.e. heating of the disk protons by the penetrating 
external protons, and internal viscous heating according 
to standard disk theory. We have neglected these here, as 
their energetic contribution is negligible in our model (the 
main energy channel is from the incident protons to the 
radiation field, via the disk electrons). If such an instabil- 
ity exists, it might heat the disk protons sufficiently for 
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the disk to expand and feed its mass into the surrounding 
ADAF. 
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